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Abstract. Grover’s algorithm is a quantum query algorithm solving the 
unstructured search problem of size N using 0{yN) queries. It provides 
a significant speed-up over any classical algorithm IGro96l . 

The running time of the algorithm, however, is very sensitive to errors in 
queries. It is known that if query may fail (report all marked elements as 
unmarked) the algorithm needs Q(N) queries to find a marked element 
IRS08| . |AB+13| have proved the same result for the model where each 
marked element has its own probability to be reported as unmarked. 

We study the behavior of Grover’s algorithm in the model where the 
search space contains both faulty and non-faulty marked elements. We 
show that in this setting it is indeed possible to find one of non-faulty 
marked items in 0(\/N) queries. 

We also analyze the limiting behavior of the algorithm for a large number 
of steps and show the existence and the structure of limiting state pum- 


1 Introduction 

Grover’s algorithm is a quantum query algorithm solving the unstructured search 
problem of size N using 0(yN) queries. It is known that any deterministic 
or randomized algorithm needs linear time (number of queries) to solve the 
above problem. Thus, Grover’s algorithm provides a significant speed-up over 
any classical algorithm. 

The running time of the algorithm (number of queries), however, is very sen¬ 
sitive to errors in queries. Regev and Schiff 1RS08 ? have shown that if query has 
a small probability of failing (reporting that none of the elements are marked), 
then quantum speed-up disappears: no quantum algorithm can be faster than 
a classical exhaustive search by more than a constant factor. Ambainis et al. 
IAB+131 have studied Grover’s algorithm in the model where each marked ele¬ 
ment has its own probability to be reported as unmarked, independent of prob¬ 
abilities of other marked elements. Similarly to the result of [RS08] . they have 
shown that if all marked elements are faulty (have non-zero probability of failure) 
then the algorithm needs f2(N) queries to find a marked element. 

* This research was supported by EU FP7 project QALGO (Dmitry Kravchenko, 
Nikolajs Nahimovs) and ERC project MQC (Alexander Rivosh) 















Although, technically the model of IAB+131 allows one non-faulty marked 
element^ (element with zero probability of failure) this case was not included 
into the analysis. 

We study the behavior of the algorithm in the model where the search space 
contains both faulty and non-faulty marked items. Specifically, we focus on the 
case where the search space contains multiple non-faulty and one faulty marked 
element. We analyze the effect of a fault on the state of the algorithm and 
show that in this setting it is indeed possible to find one of non-faulty marked 
elements in 0(y/N) queries. Up to the best our knowledge, this is the first 
demonstrgggggation of query fault modes which can be tolerated by the Grover’s 
algorithm. 

We also analyze the limiting behavior of Grover’s algorithm for a large num¬ 
ber of steps and show the existence and the structure of limiting state pu m . The 
limiting state is a mixture of a faulty marked element \i) with probability 1, the 
uniform superposition of all non-faulty marked elements with probability | and 
the uniform superposition of all non-marked elements with probability 7. Using 
the approach of IAB+131 one can show that convergence time is O(N). 

The contrast between this result and the algorithm finding a marked element 
in O(VN) steps can be explained by the probability of finding a marked element 
oscillating between 17(1) and o(l) until these oscillations finally die off after O(N) 
steps. 


2 Technical preliminaries 

We use the standard notions of quantum states, density matrices etc., as de¬ 
scribed in jKLM07j . 


Grover’s algorithm [Gro96j 

Suppose we have an unstructured search space of size N. Grover’s algorithm 
starts with a starting state \i/j 0 ) = X^iLi I*)- Each step of the algorithm 
consists of two transformations: Q and D. Here, Q is a query to a black box 
defined by 

° QV) = — |*) if * is a marked element; 

° QV) = K) if * is not a marked element. 


D is the diffusion transformation described by the following N x N matrix: 
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1 The limitation of at most one non-faulty marked element comes from the probability 
independence assumption - two or more marked elements with zero error probability 
of failure would not be independent. 













We refer to | ip t ) = (DQY \il>o) as the state of Grover’s algorithm after t time 
steps. 

Grover’s algorithm has been analyzed in detail and many facts about it are 
known |Amb04j . If there is one marked element i, the probability of finding it by 
measuring \ip t ) reaches 1 — o(l) for t = 0(\/N). If there are k marked elements, 
the probability of finding one of them by measuring \tp t ) reaches 1 — o(l) for 
t = 0(^/Njk). 


Spherical trigonometry 


Spherical trigonometry is a branch of geometry which deals with the relationships 
between trigonometric functions of the sides and angles of the spherical polygons. 
Trigonometry on a sphere differs from the traditional planar trigonometry. For 
example, all distances are measured as angular distances. 

In the context of this paper we need only a few basic formula for right spher¬ 
ical triangles. Let ABC be a right spherical triangle with a right angle C. Then 
the following following set of rules (known as Napier’s rules) applies: 


cosc 

= cos a cos b 

(m) 

tan b = 

cos A tan c 

(R6) 

sin a 

= sin A sin c 

(R2) 

tana = 

cos B tan c 

(R7) 

sin b 

= sin B sin c 

(R3) 

cos A = 

sin B cos a 

m 

tan a 

= tan A sin b 

(R4) 

cos B = 

sin A cos b 

(R9) 

tan b 

= tan B sin a 

(R5) 

cosc = 

cot A cot B 

(RIO) 



Fig. 1. Spherical trigonometry basic triangle. 


For more detailed introduction into spherical trigonometry see [Tod86| . 










3 Model and results 


Error model 

Suppose we have a search space of size N containing k marked elements i\, 12 ,..., ik- 
First k — 1 marked elements are non-faulty - the query always returns them as 
marked. Last marked element is faulty - the query might return it as unmarked. 

More formally, on each step, instead of the correct query Q , we apply a faulty 
query Q' defined as follows: 

° Q'Vk) = | ik) with probability p; 

° Q'\ik) = — |ife) with probability 1 — p; 

° Q'\j) = Q\j) if 3 7^ ik- 

Summary of results 

First we show that if there is at least one non-faulty marked element, then it is 
still possible to find a non-faulty marked element in 0(\fN) queries with (9(1) 
probability. 

Theorem 1 Let k > 3, then we can choose t = 0(i/N/k) so that, if we run 
Grover’s algorithm for t steps and measure the final state, the probability of 
finding a marked element is at least cos 2 | = 0.85.... 

For k = 2, the probability of finding a marked element is at least cos 2 | = 
0.85 ... under a promise that at most one fault occurs and at least 0.74 ... in the 
general case. 

We conjecture that, for k = 2, the probability is at least 0.85... even in the 
general case. 

Second, we analyze the limiting behavior of the algorithm for a large number 
of steps and show the existence and structure of limiting state pu m . 

Theorem 2 Let pt be the density matrix of state of Grover’s algorithm with a 
faulty oracle Q' after t queries. Then, the sequence pi , p 2 ,. ■. converges to 

Pum = ^IV’+XV’+I + \\ik)(ik\ + ^,IV ; -)(V , -I 

where \ip+) = lb) uniform superposition over all non-faulty 

marked elements and \ip~) = ^/^_ k N) uniform superposition over 

all non-marked elements. 

4 Search 

In this section we analyze the evolution of the state of Grover’s algorithm in 
presence of multiple non-faulty and one faulty marked item. First we review the 
original Grover’s algorithm, then we describe the effect of faults on the state 
of the algorithm. We derive upper bounds on the effect of faults and provide a 
modification of the original Grover’s search algorithm which finds one of non- 
faulty marked items with 0(1) probability in 0(i/~N) queries. 



4.1 No faulty marked items 


Let us first consider the very basic search problem of Grover’s algorithm. Namely, 
we have N items among which k are markec0 

Operator D is symmetric w.r.t. permutations of amplitudes of all items, and 
operator Q is symmetric w.r.t. permutations of amplitudes of marked items, 
as well as permutations of amplitudes of non-marked items. So, on any step t 
amplitudes of all marked items are equal to each other and amplitudes of all 
non-marked items are equal to each other. Thus, we can represent \ip t ) as: 

I A) = J2a t \i) + Pt\j ; ), 

i£U j£M 

where U stands for the set of indexes of non-marked items and M stands for the 
set of indexes of k marked items, at and /3* denote the amplitudes of respectively 
a non-marked item and a marked item on step t. At each step of the algorithm 
we shall take care of two numbers only: 

atV N — k and Pt.Vk. (1) 


Since | ipt) is a unit vector, we have 


ieu j€M 


= 1 and corre- 


Thus, values 0 meet the equality (a t \/N — k'j + (^fitVk) 
spond to a point on the unit circle. 

Initially all amplitudes are equal, so a = /3 = . and the numbers 0 are 


aoV N — k 


y/N — k 

Vn 


and 


/3 0 Vk 


Vk 

Vn' 


( 2 ) 


During the first step of the algorithm operator D does not change amplitudes 
of the state |^o)> an d operator Q negates amplitudes of all marked items: /?i = 
—(3 o = So the numbers ([!]) are 


oqV-N — k 


VW^k 

VN 


and 


f3\Vk 


y/k 

Vn' 


( 3 ) 


According to 0 and §, transformation |^o) IV’i) can be represented 
on the unit circle as shown on Figure [2] As before, we assume C IV, so that 
Vk <C VN — k, and the angle between | 0 O ) an d \ipi) is 


2 arcsin 




( 4 ) 


2 It is usually considered that k N, as for -A > A with sufficiently large A the 
search problem can be trivially solved by a probabilistic algorithm in time O (A^ 1 ). 
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Fig. 2. The first step of Grover’s algorithm 


(this approximation holds for small-valued -^=). 

Similarly, all further applications of operator DQ are nothing but clockwise 
rotations by angle ~ 2^|. After ~ 2v A/7]v = °’ 78 yT ^ such relations the 
resulting state 785 ^/WJk\) reaches the neighborhood of the point (0,-1). 
Measuring |^ Q 7g5 j) results in getting index of a marked item, with prob¬ 

ability almost 1. 


4.2 One Faulty Marked Item 

Let us now consider the case with 

• N — k non-marked items, 

• k — 1 marked items, and 

• 1 faulty marked item. 

To simplify the analysis we shall interpret the step of the algorithm as consequent 
application of three operators: ordinary diffusion D and ordinary query Q , and 
with probability e - error E, which negates back the amplitude of the faulty 
marked item. 

As the operation E is probabilistic one must deal not with a pure state I'i/’t), 
but with a mixed state p t (probabilistic mixture of pure states). We shall denote 
components of the mixture after t steps as | if™), where w € {0,1} ( stands for 
the sequence of t events: 0 - the query has negated the amplitude of the faulty 
marked item (DQ), and 1 - the query has left that amplitude of the faulty 
marked item unchanged ( DQE ). So the mixture pt looks as follows: 

E eH (i-e) t_H mm- 

™e{o,i }* 


Pt = 



Transformations D , Q and E are symmetric w.r.t. permutations of ampli¬ 
tudes of non-faulty marked items, as well as permutation of amplitudes of non- 
rnarked items. So, in any state | i/i™) of the mixture | ip*), amplitudes of all non- 
faulty marked items are equal to each other and amplitudes of all non-marked 
items are equal to each other. Thus, we can represent \^)f) as: 

m = E<i*> + E m+it h), (5) 

ieu j eM, 

i&k 

where U stands for the set of indexes of non-marked items and M stands for 
the set of indexes of k marked items, af, (3™ and 7 “’ denote the amplitudes of 
respectively a non-marked item, a non-faulty marked item and the faulty marked 
item. 

At each step of the algorithm for each of 2 4 scenarios w we shall take care of 
three numbers: 


a™VN - k, (3?Vk - 1 
Since | if)™) are unit vectors we have 


and 


It 


( 6 ) 


EKf + E m 2 +W ) 2 = i- 


ieu 


jeM, 

j¥=f 


Thus, values Q meet the equality (a™\/N — k) 2 + (ffiy/k — l ) 2 + (y “’) 2 = 1 
and correspond to a points on the unit sphere. 

Initially the mixture consists of state |^o) with amplitudes of all items being 
equal, so a = /? = 7 = , and the numbers ([ 6 ]) for t = 0 are 


r— - r y/N — k f- \/k — 1 

During the first step of the algorithm 


and 


7o = 


Vn' 


(7) 


D does not change amplitudes of the state |^o)i 
Q negates the amplitudes of all marked items: = 7 ° = — 

E negates back the amplitude of the faulty marked item: y-f = — 7 ° = ^=. 
So the numbers § for t = 1 are as follows: 


,;vW^E=^S, SlVTU T = -^3, 7? = -^ 

' ' Jn* ' ’I 1 i 1 ' 1 ’A ’ ~ VN 


Vn 


for w = 0 ; 
for w = 1 . 


( 8 ) 


According to ([7]) and (J 8 |, transformation |^o)(^o| D ®^ E \ (1 — e) k'i)(V , i 
can be represented on the unit sphere as shown on Figure |3| 























Fig. 3. The first step of Grover’s algorithm with a faulty marked item. Size of a 
ball corresponds to a probability of the state in the mixture. IV’oXV’ol = 1 |"0o) an d 

IV’iXV’il — (i — e ) IV’iXV’il + e l V’i > <"011 





Note that if e = 0 the state of the algorithm travels clockwise along the 
slanted ortho dromc|^] which contains points |^ 0 ) and (1,0,0). The travel lasts 
until the state reaches the neighborhood of the vertical orthodrome (hereafter 
we call it meridian) where a = 0 after t ss 0.785... yjN/k steps. 



Fig. 4. Probable routes of the initial state I'lpo) in the run of Grover’s algorithm with 
a faulty marked item. 


If k -C N the rotation angle Q is sufficiently small. The run of the algo¬ 
rithm can be viewed as clockwise rotation of a state in parallel to the slanted 
orthodrome (HQ-movement) with occasional (e-probable) up-and-down jumps 
around the horizontal orthodrome (hereafter we call it equator ), which corre¬ 
spond to operator E, as shown on Figure 0 During the first 0.785 ... y/N/k 
steps, the state cannot go out of the area which is covered with arrows on the 
figure. 

As we already mentioned, 0.785... yjN/k steps are necessary to reach the 
desired plane where a = 0, given that no fault occurs on the way (in our notation: 
in expression <§. a“-° 6 ^ « 0). 

But what could the length of the route be if some faults occur on the state’s 
way to the desired plane? In the two following subsections we will derive upper 
bounds for the effect of these faults. 


At most one fault First, let us assume that the total number of faults is at 
most one. Although this assumption seems to be rather implausible, we have 
some arguments for it: 


3 Orthodrome, also known as a great circle, of a sphere is the intersection of the sphere 
with a plane which passes through the center point of the sphere. 









o for sufficiently small e, we have ... < e 3 C e 2 C e, so that probability of 
more than one fault Y?f —2 (/) (1 — e ) 4 < t 2 e 2 could be neglected for 

number of steps t £ o Q); 

o as we shall see later, the second and all subsequent faults have smaller effect 
and even have great chances to drive the state closer to the desired meridian. 



Let us calculate the effect of a fault in the sense of its projection onto the “no¬ 
faults” route. On the Figurejfijwe illustrate e-probable transformation \ipi"°) -^> 
I'i/’j" 1 ), which happened on some step t. The fault increases the angular distance 
to the desired meridian (a = 0) by c Prr . Using rules of spherical trigonometry 
((|R8l) and (TUT])) we have: 


c err = arctan 


^tan 2 a cos arcsin 


cos A 
cos a 


= arctan 


^tan 2a 



cos 2 A 
cos 2 a 


(9) 


where A = arctan - is angle between the two equators, and a is the distance 
between and the horizontal equator. Note that a is at most A (a = A only 

when reaches the desired meridian), so 1 — — > 0. 

In equation ([ 9 ]) we assumed that \4>i"°) is located on the bottom-margin of 
the arrow-filled area of the Figure [4] (which always holds for one fault case). For 
\ipi"°) located above the bottom margin, we should calculate ([9]) for a smaller 
value of A 1 which will result in smaller value of c err . For l^” 0 ) located above the 
horizontal equator, the fault-effect c er r is negative, i.e. the resulting state | ipt" 1 ) 
is closer to the desired meridian w.r.t. the direction in parallel to the slanted 
orthodrome. Relaxing the above-mentioned assumption, we can conclude the 
following rough bound: 













Cerr < arctan tan 2 a \ 1 — 


cos 2 A 


( 10 ) 


Now, based on the inequality <[10j) we shall derive more precise bounds. 

If the number of non-faulty marked items k — 1 > 1, then the angle A = 
arctan - £ (arctan arctan (-] = (0; \]. So we have 0 < a < A < f. 

Now let us consider different values of a. If 0 < a < then is bounded 
by 


Cerr < max arctan tan 2 a \ l — 

o <a<A, \ 


cos 2 A 


7T 

~ 4’ 


(ID 


where the inequalities become equalities for A = j and a = |. 

If | < a < f, then we can follow that the state is gone far away from 

the point “ay/N — k = 1” of the unit sphere. This distance between the point 


(1,0,0) and the state | ipi"°) can be derived from the rule (R2): 


sin a sin 7 t/6 tt 

c = arcsm ——- > arcsm —-— = — 

Sind sm7r/4 4 


( 12 ) 


Since the total distance between the point “aV-N — k = 1” and any point of 
the meridian “aV N — k = 0 ” is exactly we follow that the state is at 


most ^ — j = j f ar from the desired meridian. 


From dllb and (12) we formulate the following joint conclusion: 


Corollary 1 . At least one of the following claims holds for any state \ip' t "°) with 


0 < a < A < f : 


o either the fault-effect c err < 

o or the state already at most j far from the desired meridian “a 

0 ”. 


Any number of faults Now let us use another approach to study the evolution 

of a state in the considered settings. Transformation | iff) IV’t+i) drives the 

state \iff) clockwise in parallel to the slanted orthodrome by a distance, which 
depends on the position of this state on the unit sphere. 

On Figure [ 6 ] we show a “speed” for each possible position of the state 
On the “no-faults” route (i.e. the slanted orthodrome) the speed coincides with 

that of the original Grover’s algorithm: v G = 2arcsin^(= « After 

a state jumps up, its speed decreases depending on its distance to the slanted 
orthodrome: e.g. on the parallel circle which contains point “(3y/k — 1 = 1” of 
the units sphere, its speed is v G cos A. 

We can also calculate the speed of a state w.r.t. the direction in parallel to the 
horizontal equator, i.e. the speed of its projection onto the horizontal equator v n . 

















Fig. 6. “Speed” of a state located on some meridian b. 


Obviously, among all states on some meridian b , the uppermost state 1 4>l) has the 
least speed. Distance between | ifl) and the point ll a\^N = 1” of the unit sphere, 
can be derived from angle A, distance b and rule (R6): c = arctan Distance 
between \ij^) and the slanted orthodrome can be derived from distance c, angle 
2A and rule (R2): a' = arcsin (sin 2 A sin c) = arcsin (sin 2 A sin arctan ). 
And the speed of | JpV) 


(b) 


v G cos (a) 


= V G 



sin 2 2A tan 2 b 
cos 2 A + tan 2 b 


serves as a natural lower bound for the speed of any state located on the same 
meridian b: 


v ( b) > ( b ) = v G 



sin 2 2 A tan 2 b 
cos 2 A + tan 2 b 


Projection of the speed v ( b ) onto the horizontal equator might be slightly less: 


v n ( b ) > v ( b ) cos A > 


/ sin 2 2 A tan 2 b 

1- — -, cos A. 

cos 2 A + tan" b 


(13) 


Vertical jumps leave states on the same meridians, so faults does not affect 


value b. Moving at least at speed (13), a state will pass distance § in at most 

lo 


Jo v~lb) ^ steps- From (13) we obtain an upper bound for the number of steps 
until some meridian b* (for arbitrary number of faults, i.e. for any e): 


1 


(b) 


d b< 


1 


sin 2 2 A tan 2 b 


■d b. 


1 - 


cos 2 A+tan 2 b 


cos A 


t b » < 


0 


(14) 



























We note that on the fastest “no-fault” route a state will travel exactly at 
speed v G , so that in the same many t b * steps it can reach at most (v G tb*) th 
meridian. 


From (141 we can derive the upper bound for the distance between the two 
meridians b* and v G tb »: 


rb * 


v G t b . -b* <y ^ 


1 


fl / /1 sin 2 2A tan 2 b a 


■ d& - b* 


(15) 


For example, if we fix b* = 17r, then for any 0 < A < 0.1953...7T, value 
(15) does not exceed ^. That is, when the fastest state of a quantum ensemble 


reaches meridian b* + f = |7r, the slowest state of the ensemble with certainty 
reaches at least meridian b* = |7r. 


4.3 Proof of theorem [l] 

We run standard Grover’s algorithm | times longer than usually, and then per¬ 
form a measurement. 


If at most one fault is promised, and A < 0.257t|^] then we use Corollary [I] 
and follow that any component of the resulting mixture \ip*) is at most | 
far from the meridian “a = 0”. 

If A < 0.1953 ... 7rJ^]then we substitute b* = |7r in (15) and follow exactly 


the same. 


Measurement of such \ip*) results in finding a marked item (the faulty or a 
non-faulty one) with probability at least cos 2 | = 0.853553... 

• Otherwise, if 0.1953 ... n < A < 0.257T and there is no promise on the number 
of faults, it means that there is exactly 1 non-faulty marked element (so 
A = 0.25tt). 


In this specific case we run standard Grover’s algorithm ss 1.34 times longer than 
usually, and then perform a measurement. We substitute A = 0.257T, b* ss 0.337T 
in ( fl5| and follow that value (15) does not exceed ss 0.347T. That is, when the 
fastest state of a quantum ensemble reaches meridian b* + 0.34-7T = 0.677t(= 
0.5 x 1.347r), the slowest state of the ensemble with certainty reaches at least 
meridian b* = 0.337T. Measurement of such \ip*) results in finding a marked item 
(the faulty or a non-faulty one) with probability at least cos 2 0.177T = 0.74.... 

□ 


4 A < 0.257T means that there is at least as many non-faulty marked items as faulty 
marked items. Since we limit our considerations with only one faulty marked item, 
it suffices with only one non-faulty marked item. 

5 For one faulty marked item, it means existence of at least [arccot 0.1953 ... if\ = 
[1.02047] = 2 non-faulty marked items 










5 Limiting behavior 


In this section we analyze the state of the algorithm after a large number of steps. 
We consider a density matrix of the state and describe how transformations of 
the algorithm change the state. We prove an existence of the limiting state and 
give probabilities of finding a non-faulty marked, as well as faulty marked item 
then measuring the state. 


Proof of theorem [2] 

Consider the density matrix p t of the quantum state of Grover’s algorithm after t 
queries. Due to symmetry, we can assume that the first k basis states correspond 
to the marked elements. Note that Grover’s algorithm acts in the same way on all 
non-faulty marked elements, as well as on all non-marked elements. Therefore, 
the state of the algorithm is a probabilistic mixture of pure states of the form 


a#+) + P\h) + tI'0—) - 

The density matrix p tl then, takes the form 


a 

... a 

a' c 

... c 

a 

... a 

a' c 

... c 

a' 

... a ' 

b d 1 

...d! 

c 

... c 

d! d 

... d 

c 

... c 

d! d 

... d _ 


(16) 


(17) 


because the density matrix for every pure state (161 in the mixture pt is of this 
form. 

Next we consider the effect of the faulty query Q' and the diffusion transfor¬ 
mation D on the density matrix p. 


Lemma 1 The effect of the faulty query Q' on the density matrix pt is: 


a i ^ a 

a 1 i ^ — (2p — 1 )a' 
b ^b 
c i ——c 
d' i ^ (2 p — 1 )d' 
d i—>■ d 


(18) 


Proof. Follows from IAB+131 formula (2)]. 

Lemma 2 The effect of the diffusion transformation D on the density matrix p 
is: 















Pij ^ 4U - 214,. - 2V.j + PiJ , (19) 

where V = j Pip * s tt le average of all elements of p, 14,- = ■ Pij is 

the average of i th row and V.j = 'flu Pip is the average of j th column of p. 


Proof. First we right multiply p = (pij) by 


-1 + -2 
T n 

2 

TV ‘ ‘ 

2 

2 

TV 

-1+2. 

1 ' TV ’ ’ 

2 

AT 

2 

TV 

2 

AT 

-1+ 2_ 

• 1 ' TV 


We get matrix p' = pD where 

2 

Pij = E! h >i ^ k Jf ~ P/t = 2V),. — Pip¬ 
it 


( 20 ) 


Next we left multiply (20) by D and get matrix p" 


ftp = E - Mb) - (2Vi, 


_Dp_D where 


Pij ) — 


4 _ 2 

= ]y E Vk,- ~ Jy E P fc b “ 2 ^V + P*b = 

41 / - 21 /,. - 2V.J + Pij. 


As V),. = 14,j (this follows from the structure of the density matrix p t ) we can 
use Vi to denote both 14 ,. and V.j. 

Using the same argument as in AB+131 we can prove lim^oo a! = 0 and 
lim^oo d' = 0. Consider t' such that for all t > t' inequalities a' < e and d! < e 
hold. We will prove that at this moment all elements of pt should be O(e) close 
to their final values. The proof is done in two steps. First we prove that all V) 
can differ by at most O(e). Next, as a consequence, we derive Ap i:j = 0(e). 

Claim 1 Vi,j : 14 — Vj = 0(e). 

Proof. First, note that 14 = ... = 14 _i and Vjv-fc = ... = Uv- This follows 
from the structure of p t . 

Next, consider effect of D on elements p\^ and pk,N- Both elements are O(e). 
Corresponding changes Api^ = 41/ — 214 — 214 and Apk t N = 4U — 214 — 214v 
should also be O(e). Thus, 

Ap\ } k — Apk,N = 214 — 214v = O(e) (21) 


from which follows 14 — 14v = O(e). 







Now consider Ap\^ = 4V — 2V\ — 214 = 0(e). By using (21) we can write 


V = ^((k - l)Vi + 14 + (AT - fc)Uv) = i((7V - l)Vi + V)t + O(e)). 
Putting this into Api k and opening brackets we get 


27V — 4 27V — 4 ^ 

-Pi- 77 —Vk = O(e) 


TV 


TV 


from which follows Vi — 14 = O(e). 

Corollary 2. D changes elements of pt by at most O(e). 

Proof. The change of element p l :) is equal to Ap l:l = 4V — 2V) — 2V). Using 
claim [l] we can rewrite this equation as 

Apij = — (TVV) + O(e)) — 4V) + O(e) = O(e). 


Corollary 3. c = O(e). 


Proof. Consider the effect of Q' on V\ = (k — l)a + (TV — fc)c + O(e). Q' leaves 
a unchanged, but changes the sign of c. As the claim [T] should hold for all 
subsequent steps, Q' can not change value of Vi by more than 0(e). Thus, 
Vi(c) — Vi(—c) should be O(e), which means that c = O(e). 


Corollary 4. Q' changes elements of p t by at most O(e). 


Proof. Trivially follows from (18) and the previous corollary. 


We can chose e to be arbitrary small, thus, we have 


Plim — 


a 

... a 0 0 ... 

O' 

a 

... a 0 0 ... 

0 

0 

... 0 bO ... 

0 

0 

...00 d... 

d 

_0 

... 00 d... 

d 


As ( k — l)a = b = (TV — k)d (follows from Vi = 14 = Uv) and {k — l)a + b + 
(TV — k)d = 1 (property of density matrix) we have ( k — l)a = b = (TV — k)d = | 
which completes the proof of the theorem. 

□ 
















6 Summary and open problems 


In this paper we focus on the case where search space contains multiple non- 
faulty and one faulty marked element. 

First we show that if there are at least two non-faulty marked elements or 
it there is at most one fault, then it is still possible to find a non-faulty marked 
element in 0(yN) queries with ©(1) probability. 

Second, we analyze the limiting behavior of the algorithm for a large number 
of steps and show the existence and structure of limiting state pu m . 

It is an open question to generalize the results for more than one faulty 
marked element. Although the generalization of Theorem [2] seems to be straight¬ 
forward it is not given here. The generalization of Theorem [l] might be tricky as 
one needs to deal with hyper-spherical geometry. One simple case, in which the 
theorem still holds for multiple faulty marked elements, is synchronized-fault- 
case, where each query either fails or succeeds on all faulty marked elements. 
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